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Abstract 

As photodisintegration of deuterons mimics the disintegration of deuterons by neutrinos, the accurate measurement of the 
radioactivity from thorium and uranium decay chains in the heavy water in the Sudbury Neutrino Observatory (SNO) is 
essential for the determination of the total solar neutrino flux. A radium assay technique of the required sensitivity is described 
that uses hydrous titanium oxide (HTiO) adsorbent on a filtration membrane together with a j8 — a delayed coincidence 
counting system. For a 200 tonne assay the detection limit for 232 Th is a concentration of ~ 3 x 10~ 16 g Th/g water and for 
238 U of ~ 3 x 10~ 16 g U/g water. Results of assays of both the heavy and light water carried out during the first two years 
of data collection of SNO are presented. 
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1. Introduction 



The Sudbury Neutrino Observatory (SNO) [1] is a 
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at a depth of 6800 ft in INCO's Creighton mine at 
Sudbury, Canada. Its main aim is to observe neutrinos 
from the sun and investigate the origin of the observed 
deficit in their flux (the Solar Neutrino Problem). The 
detector uses ^9500 photomultiplier tubes (PMTs) on 
a geodesic sphere of diameter 17.8 m to observe the 
Cherenkov light produced as a result of neutrino in- 
teractions occurring in the 1000 tonnes of D 2 0. The 
D 2 is held in a 12 m diameter acrylic sphere, which 
is surrounded by a shield of 7000 tonnes of H z O con- 
tained in a 34 m high barrel-shaped cavity of maxi- 
mum diameter 22 m. 

SNO can measure the flux of electron neutrinos 
from the sun via the charged current (CC) reaction: v e 
+ d^p + p + e~, and the total flux of all active neu- 
trino flavours via the neutral current (NC) reaction v x 
+ d — > p + n + v x , (x = e, jj,, t). The first results from 
SNO on the total flux of all active neutrinos have re- 
cently been published [2], in which a comparison of 
the CC and NC flux provides evidence that neutrino 
flavour changing is the cause of the Solar Neutrino 
Problem. 



2. Radioactivity limits in SNO 

Of the naturally occurring radioactive nuclides, it is 
only in the decay chains of the two heavy nuclides, 
232 Th and 238 U, that gammas with energies greater 
than 2.22 MeV can be produced. These can photodis- 
integrate the deuteron and hence mimic a NC event, 
as the signature of the NC reaction is the production 
of a free neutron. (At 6800ft depth, the number of 
cosmogenically-produced radioactive nuclides is neg- 
ligible.) 

The number of neutrons predicted by the Standard 
Solar Model that are produced by the neutral current 
interaction of solar neutrinos with deuterons is about 
0.015 per day per tonne of heavy water. In the tho- 
rium decay chain 208 T1 is produced with a branching 
ratio of 36% and in its decay a 2.614 MeV gamma is 
emitted with a 99.8% probability. The chance that this 
gamma photodisintegrates a deuteron is 2.07xl0~ 3 , 
giving a probability that a neutron is generated of 
7.45 xlO -4 per Th decay. An upper limit on 232 Th 
of 3.8xl0 _15 g/g (throughout this paper, g/g means 
gram of 232 Th or 238 U per gram of D 2 or H 2 0) cor- 



responds to a photodisintegration background of one 
neutron per day per 1000 tonnes of heavy water. This 
level is less than 10% of the predicted NC rate and 
corresponds to an average Th decay rate of about one 
decay per day per tonne of heavy water. 

In the uranium decay chain 214 Bi is produced, and in 
its decay the probability that a gamma with an energy 
greater than 2.22 MeV is emitted is only ~2%. The 
probability that a neutron is generated is 3.12xl0~ 5 
per U decay and this sets a corresponding upper limit 
on 238 U of 3.0x 10~ 14 g/g. The light water surrounding 
the heavy water must also be very pure to reduce ra- 
dioactivity backgrounds though the requirements are 
not so severe, with the concentration limits for 232 Th 
and 238 U, based on early simulations of the detector, 
being 37xl0~ 15 g/g and 45xl0~ 14 g/g, respectively. 

To achieve this level of radiopurity, both the heavy 
and light water in SNO were purified via several stages 
including filtration, ion exchange and reverse osmo- 
sis, which are described in detail in [1]. The water can 
be recirculated for further purification and assay. All 
the materials used in the construction were also se- 
lected for low radioactivity. To reduce any radon gas 
contamination of both the light and heavy water, all 
components in the water systems were chosen for low 
radon diffusion and emanation and the surfaces of the 
D z O and H 2 are covered with N 2 gas obtained from 
liquid nitrogen. 

It is essential to be able to assay the heavy water 
in the SNO detector at these very low levels of ra- 
dioactivity to check that the required radiopurity has 
been achieved and measure the NC background. With 
this aim in mind, the hydrous titanium oxide (HTiO) 
ion-exchange system was developed to assay Ra, Th 
and Pb isotopes. It is one of four complimentary tech- 
niques used in SNO: the other three being an MnO x 
assay method, for Ra isotopes, degassing for 222 Rn 
from the U chain, and direct counting of the amount 
of Cherenkov light from the decay of 208 T1 and 214 Bi. 
These other techniques will be described elsewhere. 



3. The HTiO assay method 

In this technique heavy or light water is passed 
through HTiO trapped on the filtration fibres of mi- 
crofiltration membranes with a pore size of 0.1 jUm. 
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(Hydrous titanium oxide is a non-stoichiometric com- 
pound so the chemical formula HTiO is only repre- 
sentational.) Hydrous titanium oxide is an inorganic 
ion-exchanger [3,4] that has been used in the nuclear 
industry for the extraction of heavy ions. It is formed 
as a white colloidal suspension by the hydrolysis of 
titanium salts [5], which forms particles of sizes rang- 
ing from 1 to 100 microns, with a mode at around 10 
microns. Examples of its use are given in [6,7,8]. The 
initial investigations for SNO on the use of HTiO for 
extracting radium, thorium and lead are described in 
[9,10,1 1]. In this paper the extraction and assay of ra- 
dium from the Th and U chains is reported. 

In the HTiO assay procedure there are five steps: 
deposition, extraction, elution, secondary concentra- 
tion and counting. The HTiO ion-exchanger is first de- 
posited onto microfiltration membranes: 0.25 m long 
filters for a light water assay and 1.0 m long filters 
for a heavy water assay. Then the columns contain- 
ing the loaded filters are taken underground to ex- 
tract 224 Ra from the Th chain and 226 Ra from the U 
chain from a large volume of D 2 (^200 tonnes) or 
H 2 (^30 tonnes). After extraction, the columns are 
brought back up to the surface laboratory, and the ra- 
dium is eluted with 15 I of dilute nitric acid. In the 
secondary concentration stage, the radium in the 15 / 
nitric eluate is further concentrated down to ^10 ml 
of final eluate. This is then mixed with ^40 g of liq- 
uid scintillator and counted using a j3 — a delayed- 
coincidence liquid-scintillation counter. 

The measurement of the concentration of 224 Ra and 
226 Ra therefore requires the determination of the ex- 
traction efficiency of the loaded filters for a through- 
put of ^200 tonnes of water, the subsequent elution 
and secondary concentration efficiencies, and finally 
the detection efficiency of the j3 — a counters. The in- 
ference of the concentration of the elements in the Th 
and U chain in the heavy or light water then requires 
knowledge of the equilibrium conditions within the 
chains, and of the extent of any radioactive species 
plating out within the water systems. 

Initially only heavy water assays were made, us- 
ing hollow fibre 1.0 m filters, and then the assays 
were extended to both the light and heavy water using 
0.25 m and 1.0 m pleated membrane filters, respec- 
tively. In this paper the determination of the extrac- 
tion and elution efficiencies for both these filters are 
first explained; then the j3 - a counting system and 



its detection efficiency are described; finally the assay 
results are presented and discussed. 

3.1. Production of HTiO 

The HTiO is produced by adding 400 ml of 15% 
w/v Ti(S0 4 ) 2 to 3 I of ultrapure water (UPW) and then 
adjusting to pH 12 with 3.5 I of 0.5 M NaOH. Af- 
ter allowing the solution to sit in a cool dark area for 
one week, it is rinsed with UPW to reduce the Na + 
and S0 4 2 ~ contaminants. The first step in rinsing the 
HTiO involves reducing the volume of the undisturbed 
HTiO solution to ^5 I by aspirating the supernatant. 
The remaining solution is then stirred for 5 minutes 
and equally divided into eight centrifuge bottles and 
centrifuged at 3400 rpm for 7 minutes. After aspirat- 
ing and discarding the supernatant, the precipitate is 
combined with ^600 ml of UPW and stirred for 5 min- 
utes. The centrifuging and rinsing steps are repeated 
two more times. In the final step, the precipitates in 
the eight centrifuge bottles are combined in an amber 
bottle and the weight of the HTiO solution is made up 
to 5 kg with UPW 

3.2. Hollow fibre filter system 

Before 2001, only heavy water assays were made 
using two hollow fibre 1.0 m Amicon H26 filters in 
parallel, each with a pore size of 0.1 /im and a sur- 
face area of 2.4 m 2 . About 1.25 g Ti in the form of 
HTiO was deposited onto each filter, and the coverage 
was therefore 0.5 g Ti/m 2 . A specially constructed de- 
position and elution rig in the surface laboratory was 
used for this purpose. Because of the need to keep ra- 
dioactivity backgrounds to a minimum this rig is made 
principally out of low activity polypropylene. 

The HTiO is deposited by adding the required 
amount of concentrated HTiO stock to a reservoir on 
the elution rig. The volume is then made up to 15 / 
with UPW and circulated through the filter at ^85 
//min for 10 minutes. 

To determine the radium background from the 
HTiO, the chemicals used in the elution and sec- 
ondary concentration, the filters and the elution rig, 
the loaded filters are first eluted, without any extrac- 
tion of radium from the heavy water, and the eluate 
concentrated and counted. The same filters are then 
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Table 1 

Summary of efficiencies for the assay of radium 



Efficiency Hollow Pleated 



Extraction e ext 


95±5% 


95+5 %B 

— 10 


Elution e du 


70±5% 


80±10% 


Secondary cone e conc 


43±7% 


71±10% 


HI extraction 


95±5% 




HI elution 


53±5% 




MediaKap-10 extraction 


95±5% 


95±5% 


MediaKap-10 elution 


90±10% 


75±10% 


Chemical e ex ,-e elu -e conc 


29±5% 


48±11% 


Counting e coum 






Th chain 


45±5% 


45±5% 


U chain 


60±10% 


60±10% 


Total £exf£ e i u -£conc' £ count 






Th chain 


13±2% 


24±5% 


U chain 


17±4% 


32±8% 



a This is for heavy water; extraction efficiency 
for light water assays is 96±4%, see text. 



cleaned by dissolving the HTiO in 0.5 M HC1, by 
circulating the acid for forty minutes and leaving the 
filters to soak overnight. This is followed by a UPW 
rinse before reloading the filters with HTiO from the 
same batch. The filters are then taken underground. 

Once underground the two filters must be deuterated 
to maintain the isotopic purity of the heavy water. The 
two filters are first slowly filled from the bottom with 
~ 10 kg of D 2 over 2 hours, displacing any H 2 out 
of the top. The filters are left overnight to allow time 
for H 2 deep within the filter to diffuse out into the 
surrounding D 2 and the flow is resumed the next 
morning. After —9 hours flow and —20 kg more D 2 0, 
the density of the water coming out of the top of the 
filters is 1.102-1.103, and the deuteration of the filters 
is then stopped (the density of pure D z O is 1.105 g/cc 
at 20 °C). 

3.2.1. Extraction efficiency 

The extraction efficiency for a loaded filter depends 
on the volume of water passed per unit area, the purity 
of the water, the flow rate per unit area, the coverage 
of the adsorbent (g/m 2 ) and on the distribution or par- 
tition coefficient of the adsorbent (K d ). Provided the 
volume of water that has passed through the filter is 
such that the ratio of the amount of activity on the ad- 
sorbent to its concentration in the liquid is much less 
than K A , then the extraction efficiency will be close 



to 100%. Other contaminants in the water, if of suf- 
ficient concentration, can reduce the efficiency of the 
adsorbent if large volumes of water are processed. 

Measurements of the extraction efficiency of the 
hollow fibre filters were carried out at Oxford [10] 
using radium spikes. A small Amicon HI filter was 
used of area 0.015 m 2 and with the same pore size 
and with the same coverage of HTiO as for the large 
1.0 m filters. Using the same flow rate per unit area, 
these tests indicated that the extraction efficiency was 
—95%, even for an equivalent processed volume of 
1 kilotonne (= 200 tonnes/m 2 )and flow rate per unit 
area of —40 Z/min/m 2 . Typical values for the heavy 
water assays with the hollow fibre filters were 200 
tonnes (= 40 tonnes/m 2 ) at 20 Z/min/m 2 . At the site of 
the experiment, stable barium, for which HTiO has a 
similar but slightly reduced affinity [13] compared to 
radium, was added to check the extraction efficiency 
after an assay and the barium extraction efficiency was 
found to be -90%. 

These measurements give an extraction efficiency 
with the 1.0 m hollow fibre filters of 95±5% for Ra 
using the Ba data as the lower bound and 100% as the 
upper bound. 

3.2.2. Elution and secondary concentration 
efficiencies 

The elution efficiencies for the heavy water assays 
using the hollow fibre filters are based on 226 Ra spike 
experiments that were carried out on 0.25 m Memtrex 
pleated membrane filters at Laurentian University. It 
is assumed that the elution efficiency depends on the 
acid strength and contact time and further, because the 
efficiency depends on the binding ability of Ra with 
HTiO, that there is no significant difference in the elu- 
tion efficiencies for the hollow fibre and pleated mem- 
brane filters with the same contact time. A 20 minute 
circulation of 15 I of 0.03 M HN0 3 was used for the 
hollow fibre filters, which achieved an elution effi- 
ciency of 70%. For comparison, the elution efficiency 
of Ba, which would be slightly higher than that of Ra 
due to the relatively weaker affinity of Ba to HTiO, 
was found to be 75±10%. 

The secondary concentration of the 15 I nitric elu- 
ate was achieved by first adding —85 mg of Ti in the 
form of Ti(S0 4 ) 2 solution, adjusting the pH to 9 by 
titration with NaOH to form HTiO, with which the 



5 



Ra co-precipitates. The 15 I solution was then filtered 
through a small Amicon HI filter to trap the precip- 
itated HTiO+Ra, and the Ra eluted off the HI filter 
by using 140 ml of 0.03M HN0 3 twice for 20 min- 
utes each. About 50 mg of Ti in the form of Ti(S0 4 ) 2 
was added to the combined HN0 3 eluate. The pH was 
adjusted to 9 by titration with NaOH to form HTiO 
with which the Ra co-precipitates. Then the titrated 
eluate was passed through a MediaKap-10 filter to ex- 
tract the Ra. Finally the Ra was concentrated down 
to ~10 ml of acid solution by eluting the radium off 
the MediaKap-10 filter with 5 ml of 0.5 M HC1 used 
twice. 

The efficiency of this secondary concentration was 
determined by 224 Ra and 226 Ra spike experiments, 
carried out on HI and MediaKap-10 filters under the 
same conditions as in the actual assay. The final figure 
of 43 ±7% is the product of the extraction and elution 
efficiencies of the HI filter and MediaKap-10 filters 
that are used in this stage of the chemical procedure 
(see Table 1 for details). 

3.3. Pleated membrane filter 

From 2001, two 1.0 m Memtrex pleated membrane 
filters [14] (Osmonics, Inc., USA) in parallel were 
used instead of the hollow fibre filters as the pleated 
membrane filters had a higher pressure rating, which 
would allow higher flow rates to be used. The pleated 
membrane filter has 0.52 m 2 of surface area for a 0.25 
m filter, and therefore 2.08 m 2 for a 1 .0 m filter, with a 
0. 1 jxm pore size. The filters are cylindrical (diameter 
2.7") with a double O-ring seal at one end, and are held 
in a cylindrical holder (the column). The water enters 
at one end of the column, flows around the outside and 
through the HTiO loaded membrane into the inside of 
the filter, and out the other end of the column. 

Several modifications, mainly in the secondary con- 
centration stage where the HI filter is replaced with 
three MediaKap-10 filters in parallel, were made to 
improve the chemical procedural efficiency. 

3.3.1. Extraction efficiency 

To examine the full-scale extraction efficiency, two 
kinds of experiments were conducted: one using the 
underground water systems and the other, the facility 
at Carleton University where the extraction of radioac- 



tive spikes at the same flow rate per unit area used 
underground but with a lower mass/m 2 of ~ 1 tonne 
water/m 2 could be measured. 

Two light water assays, LWA2 and LWA3, were 
first carried out to examine the extraction in the actual 
underground light water system. In these experiments 
two 0.25 m filters with HTiO coverage of0.5gTi/m 2 
were placed underground in series to process about 
47 tonnes (=90 T/m 2 ) of light water at the flow rate 
of about 18 Z/min (=35 Z/min/m 2 ). The radioactivity 
extracted by these two columns was then separately 
measured. If the extraction efficiency is high, there 
should be little radioactivity on the downstream col- 
umn. However, in both experiments, significant activ- 
ity was found to penetrate the upstream filter and reach 
the downstream one (see Table 2). It was suspected 
that the HTiO coverage on the pleated membrane fil- 
ter might be too low, either to provide complete cov- 
erage or to cope with contaminants in the light water 
system, or that the high water flow rate might result 
in incomplete extraction. 

To solve this problem, a heavier HTiO loading (1 .25 
g Ti in the form of HTiO, equivalent to 2.5 g Ti/m 2 ) 
and a slower flow rate ~10 Z/min (=19 Z/min/m 2 ) were 
used in the light water assays LWA4-8. As a result, 
the average radium extraction of the upstream filter 
was improved to 80 ± 4%, and for the largest tonnes 
water/m 2 processed, LWA8, the efficiency was 80 ± 
5%. To improve the signal in the light water assays, the 
activity on both the upstream and downstream filters 
was added together, which gives an overall extraction 
efficiency of 96 ± 4%, taking the efficiency of both 
filters as 80%. 

After an extraction a loaded filter was examined to 
check the HTiO coverage. The passage of a large vol- 
ume of water had not caused the coverage of HTiO to 
become non-uniform and the extraction efficiency of 
80% for the upstream pleated filter, rather than 95% 
as found for the hollow fibre filter, is probably caused 
by some suspended non-electrolytes in the light wa- 
ter. These would cover the exchange sites of the HTiO 
adsorbent and therefore reduce the exchange capacity 
of the HTiO. The level of such suspended particles in 
the heavy water is expected to be much lower, because 
the heavy water is in contact with far less and cleaner 
material than the light water. (The level of any inor- 
ganic contaminant in the light water is less than 1 ppb 
and the level of organics in H z O is also low and less 
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Table 2 

The extraction efficiency for pleated membrane filters in the underground water systems 

Expt T/m 2 //min g Ti Upstream filter" Downstream filter Extraction 

/m 2 /m 2 (10" 15 g/g) (l(T 15 g/g) Efficiency* (%) 











Th 


U 


Th 


U 


224 Ra 


226 Ra 


LWA2 


88 


35 


0.5 


19.2±2.7 


0.94 ±0.21 


17.5 ±2.4 


0.99 ±0.26 






LWA3 


90 


35 


0.5 


22.0±2.6 


1.87 ±0.76 


11.8 ±2.5 


0.79 ±0.74 






LWA4 
LWA5 


21 
36 


19 
19 


2.5 
2.5 


22.1 ± 1.1 
60.1 ±5.4 


8.5±5.2 
3.73 ±0.90 


+ 6.2 
u -0 

18.6±4.2 


2 Q + 4.6 
1.11 ±0.79 


100+° 28 
69±8 


76 + 24 
' -56 

70 ±22 


LWA6 
LWA7 


66 
53 


19 
19 


2.5 
2.5 


47.7 ±3.8 
71.3±3.9 


4.37 ±1.25 
2.79 ±0.42 


10.4±3.0 
13.5±3.0 


UD1 -0.61 
O1+0.34 
UJJ -0.33 


78±6 
81±4 


86 + 14 

° -20 
88± 12 


LWA8 


82 


19 


2.5 


47.4 ±2.9 


2.22 ±0.47 


9.7 ±2.5 


43+0.48 
u -'"_0.43 


80 ±5 


81 ± 19 


RDSA11 


64 


19 


2.5 


1.07 ±0.14 


0.28 ±0.09 


+ 0.17 


01 + 009 
uul -0.01 




96 + 4 



" The equivalent Th(U) concentration at the end of the extraction. 
b Extraction efficiency = 1 - Downstream Th(U)/Upstream Th(U). 



than 30 ppb. These levels are unlikely to reduce the 
extraction ability of the HTiO for any activity.) 

To measure the extraction efficiency for heavy wa- 
ter, an assay (RDSA11) of the salinated heavy wa- 
ter (0.2% NaCl) was made using two pairs of 1.0 m 
pleated membrane filters in series. To determine the 
radium background, each pair was first loaded with 
HTiO which was then eluted, the eluate concentrated 
and the background radium activity counted. The fil- 
ters were then cleaned by circulating 0.5 M HC1, 
reloaded with HTiO and taken underground for the 
assay. The results for this assay are given in Table 2 
and show that the extraction efficiency of the pleated 
membrane filter is 100 + j 6 % for 64 tonnes of water/m 2 
and a flow rate of 19 Z/min/m 2 . 

In this experiment the 224 Ra background of the up- 
stream pair was anomalously high, indicating some 
contamination (probably from the elution rig). The 
cleaning of the filter would have cleaned the elution 
rig as well and so the average rate from four back- 
ground runs has been assumed for the upstream value 
of the 224 Ra background. The extraction efficiency is 
in agreement with 224 Ra spike experiments that were 
conducted at Carleton University to check the extrac- 
tion efficiency of radium in 0.5% NaCl solution using 
a single 0.25 m filter with a loading of 2.5g Ti/m 2 . 
Satisfactory extraction efficiencies (94%, 97%, 94%) 
were obtained at a flow rate/m 2 of 23 Z/min/m 2 with 
water amounts of 0.35, 0.58 and 0.92 tonne/m 2 , re- 
spectively. 

Combining these results, the extraction efficiency 
of the pleated membrane filter is taken to be 95 + ? % 



for all assays of SNO heavy water, whether salinated 
or not. 

3.3.2. Elution and secondary concentration 
efficiencies 

After an assay the columns are transported to the 
surface laboratory and mounted on the elution rig. To 
improve the elution efficiency a 40 minute circula- 
tion of 15 I of 0.03 M HN0 3 is used, which achieves 
an elution efficiency of 80%, compared to the 70% 
achieved previously for the hollow fibre filters with a 
20 minute circulation. 

The secondary concentration is carried out by first 
adding Ti(S0 4 ) 2 and then NaOH, as with the hol- 
low fibre filters, to form HTiO with which the Ra co- 
precipitates. Then, in order to improve the secondary 
concentration efficiency, the HI filter followed by a 
MediaKap-10 filter, used with the hollow fibre filters, 
is replaced with three MediaKap-10 filters in parallel 
through which the 15 I solution is then filtered. The 
radium is eluted from each with 3.5 ml 0.2 M HN0 3 
for 15 minutes. These are combined to form the first 
^10.5 ml eluate. The MediaKap-10 filters are eluted 
twice more to ensure all the activity is collected, and 
combined to form the second and third ~40.5 ml dil- 
ates. 

The elution and secondary concentration efficien- 
cies were obtained through 226 Ra and 228 Th spike ex- 
periments at Laurentian University. For the elution the 
same acid strength and contact time were used as in 
an assay, and for the secondary concentration the con- 
ditions were the same as in actual assays. 
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The elution and secondary concentration effi- 
ciencies are estimated to be, respectively, 80±10% 
and 71 ±10% for radium. (MediaKap-10 extraction 
efficiency is 95 ±5% and the elution efficiency is 
75±10%, see Table 1). The elution efficiencies for 
228 Th and 212 Pb were determined by 228 Th spike ex- 
periments to be very low (1.6% and 5%, respectively). 

232 Tha 1.4xl0 ,0 y 
228 Ra,3 5.8y 

228 Ac ,3 6.2/z 

228 Tha l,9y 
224 Rara 3Jd 
220 Rn« 56* 
216 Po a 0.1 5s 
212 Pb jS I0.6h 

212 Bi « 36% (S 64% 61m 
208 T1 (i 3m 212 Po a 0.3/js 

208p b 

238 U a 4.5x1 9 y 
234 Th/3 24d 

234 Pa (j 1. 2m 

234 U a 2.5xl0 5 ;; 
230 Tha 7.5xl0 4 y 
226 Ra a 1600^ 
222 Rna 3.8 d 
218 Po a 3 m 
214 Pb (3 21m 

214 Bi a 0.02% (3 >99% 20m 
210 T1 (S 1.3m 214 Po a 164/j.i 
210 Pb p 22y 

210 Bi /S 5rf 

2l0 Po a I38d 

206p b 

Fig. 1. Thorium and Uranium decay schemes 



3.4. fi — a counter and counting efficiency 

For a 200 tonne assay of heavy water with con- 
centrations of 232 Th of 3.8 x 10~ 15 g/g and of 238 U of 
3 .Ox 10~ 14 g/g, corresponding to background neutron 
production rates from deuteron photodistintegration of 
~10% of that predicted by the Standard Solar Model, 
the radium decay rates, if 50% is extracted, eluted 
and concentrated, are ~130 decays per hour (dph) for 
226 Ra , while only ^6 dph for 224 Ra. This very low 
level of activity requires counters with high efficiency 
and low background. To meet these requirements, a 
P — a counting system was developed in Oxford [10], 
which uses delayed coincidence between two sequen- 
tial decays, which follow the decay of Ra and occur 



towards the end of the Th and U chains, to reduce the 
background, together with pulse shape discrimination 
[16]. The decay schemes are shown in Fig. 1. 

For the Th chain, the j3-decay of 212 Bi and the a- 
decay of the 212 Po are utilized, while for the U chain it 
is the j3-decay of 214 Bi and the a -decay of 214 Po. The 
different half-lives of the two Po isotopes (300 ns for 
212 Po and 164 /is for 214 Po) allow these j3 — a sequen- 
tial decays to be easily identified and separated from 
each other. In previous j3 — a coincidence techniques 
[17,18], the activity was deposited as a thin film. How- 
ever, such techniques have quite low efficiencies and 
producing thin films from HTiO eluates was found to 
be difficult. To overcome these limitations, the activ- 
ity is concentrated down to a small volume of aqueous 
solution, which is then mixed into a liquid scintillator 
and the j3 — a decays counted with a photomultiplier. 

Plastic jars with their high radiopurity and low cost 
proved the ideal disposable container for the j3 — a 
counting. After testing various combinations [10], the 
final procedure chosen is that a ~10.5 ml aqueous 
sample containing the eluted radium is mixed with 42 
g of Optiphase Hi-Safe 3 liquid scintillator and placed 
in a 60 ml polymethylpentene jar. The 5 cm diameter 
jar and a 5 cm photomultiplier (an Electron Tubes 
Ltd 9266XB PMT) are optically coupled with silicone 
grease to maximise light collection. The counters are 
shielded against soft room gammas by a 2.5 cm thick 
oxygen-free, high conductivity copper housing. The 
system currently uses eight of these counters. 

The electronics block diagram is shown in Fig. 2. 
All events with coincidences within a time window 
<700 jUs are captured by a CAMAC system and their 
TDC and 3 ADC values (beta full-pulse, alpha full- 
pulse and alpha tail-pulse) are transferred to a PC. The 
ADCs are operated as QDCs and the delay for the 
alpha tail is set so that the integrated charge in the tail 
of the pulse is approximately a third of that in the full 
alpha pulse. 

A time cut on the time T between the j3 and the 
a signal is first used to tag an event as thorium-like 
(T <1500 ns, corresponding to the 300 ns half-life of 
212 Po) or uranium-like (10 fis< T <700 /is, corre- 
sponding to the 164 /is half-life of 214 Po). In order to 
improve the background rejection from random coin- 
cidences, three off-line software cuts are applied, il- 
lustrated in Fig. 3. A j3 energy cut, which is normally 
fixed at 75 keV < /3 < 2500 keV, removes cross-over 
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Fig. 2. Electronics block diagram of the j8 — a counters. 



coincidences between different counters. Alpha en- 
ergy (6 MeV < a < 14 MeV) and PSD (pulse shape 
discrimination, which is the ratio of a tail charge/ full 
charge) cuts greatly help to reduce the background 
events from random y— y coincidences. These cuts are 
set using spiked calibration sources. Different mea- 
surement conditions, such as temperature and cloudi- 
ness of sample, can lead to a slight a peak shift. There- 
fore, a and PSD cuts are adjusted by software for ev- 
ery individual sample. In Fig. 3, it can be seen that 
the true and background events can be easily distin- 
guished. 

A good linear relationship between counting rate 
and sample activity was found using blank and spiked 
sources of order 0. 1 mBq to 1 Bq of 228 Th and 226 Ra, 
and the blank source (a mixture of 10 ml 0.5 M HC1 
and 42 g scintillator) gave ^0.03 cph (counts per hour) 
for the Th chain and ~0.3 cph for the U chain ac- 
tivities. These count rates are much lower than the 
the total radium backgrounds using the pleated mem- 
brane filters for both the Th and U chains of typically 
0.50 and 1.34 cph, respectively, which correspond to 
- 8 x 10~ 16 g Th/g and to - 5 x 10~ 16 g U/g for a 200 
tonne assay. 

Efficiency calibrations were also made. The result 
of many calibrations of similar sets of counters gave 
counting efficiencies of 45±5% for the Th chain and 
60±10% for the U chain. In addition, as 228 Th (-1 
Bq) and 226 Ra (~0.5 Bq) sources are regularly used to 
check the performances of the counters, any change in 




500 1000 1500 

a channel 



Fig. 3. Thorium chain j8 — a coincidence data for assay RDA5. 
a: The full charge spectrum for the alpha (second) pulse of the 
coincidence, b: The spectrum for the PSD ratio of the tail to 
the full charge of the alpha, c: Two dimensional PSD— vs— alpha 
spectrum. The j3 — a events pass cuts on the full alpha charge and 
the PSD ratio and are plotted darker than the background events 
in all three spectra. 

the counting efficiencies can easily be seen. In some 
of our heavy water assays, a few samples had a light 
yellow colour that might be attributed to manganese 
and organic matter in the water. To check for any loss 
in efficiency, small 212 Pb spikes were added to the 
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yellow samples after counting; no significant loss in 
counting efficiency was found. 



4. Description of assay experiments 

The experiments included in this paper are named 
according to a code system which is: RDA# (Radium 
D 2 Assay), RDSA# (Radium D 2 Salt Assay- an 
assay of the D 2 with 0.2% NaCl added), LWA# (Ra- 
dium H 2 assay). From October 1999 to November 
2001, nine heavy water and eight light water assays 
have been conducted. Two 1 .0 m filters in parallel were 
used for all heavy water assays, two 0.25 m filters in 
series for all light water assays. The hollow fibre fil- 
ters were used for RDA2, RDA3, RDA4 and RDA5. 
For the other assays, the pleated membrane filters were 
used. 

The accurate measurement of the amount of radium 
in the water volume also requires the determination 
of any contribution from (a) the pipes connecting the 
acrylic vessel and cavity to the filters, and from (b) the 
filters, the elution rig, the secondary concentration ap- 
paratus, the chemicals and the counters. To determine 
(a) an assay was carried out of water which was circu- 
lated through the flow and return pipes to the acrylic 
vessel but with the vessel by-passed. The result was a 
limit for Th of less than 0.6x 10~ 15 g/g [12]. To deter- 
mine (b) a background measurement using the same 
filter is always carried out a few days before an assay 
with the same procedure but no extraction step. 

To minimise the procedural background, the filters, 
the counting jars and the tubing used in the secondary 
concentration are only used once for each assay. The 
titanium concentration in all solutions produced in the 
entire procedure was measured by inductively coupled 
plasma mass spectrometry and UV-visible spectrom- 
etry to check for any anomalous Ti loss; none was 
found. Since May 2001 (assay LWA4 and subsequent 
assays), the HTiO coverage on the pleated membrane 
filters has been increased from 0.5 to 2.2-2.5 g Ti/m 2 . 

The heavy water assays sampled the D 2 by draw- 
ing water from the bottom of the acrylic vessel and 
returning it to the top. The light water assays sampled 
the H 2 between the acrylic vessel and the photomul- 
tipliers, where the radioactivity contributes most to the 
background of the neutrino-induced signal. 



4. 1 . Data Reduction 

In the Th chain (Fig. 1), the rate of delayed coinci- 
dences is directly proportional to the specific activity 
of 212 Bi, X Bi N Bi , where X is the decay constant, and 
N is the number of atoms. If the radionuclides above 
212 Bi are not in secular equilibrium, then the rate of 
coincidences varies with time as the chain equilibrates. 
A spectrum of the counting rate plotted against time 
can be fitted to competing exponentials, and the spe- 
cific activity of isotopes further up the chain can be 
derived. In this fashion the activity levels of 228 Th, 
224 Ra, 212 Pb and 212 Bi can be measured. Likewise, in 
the U chain (Fig. 1), the rate of 214 Bi decays, plot- 
ted against time, can be used to fit the activity levels 
of 226 Ra, 222 Rn, 214 Pb and 214 Bi. The fitting program 
uses a log likelihood method [10] with Poisson dis- 
tributed bin heights. 

An example of spectrum fitting is shown in Fig. 4. 
The total j3 — a coincidence rate is given with the con- 
tributions from different radionuclides in the chain. 
For the Th chain data in Fig. 4, the curve labelled 
224 Ra shows the contribution from the activity of 224 Ra 
alone; initially the j3 — a coincidence rate builds up 
in a time determined by the 212 Pb lifetime (10.6 hr) to 
a maximum 37 hr after the start of counting. It then 
decays away at a rate determined by the half life of 
224 Ra (3.66 d). For the U chain data shown in Fig. 4, 
it should be noted that 214 Pb and 214 Bi can not be 
separated very well due to their close half-lives (26.8 
min for 214 Pb and 19.9 min for 214 Bi), and their con- 
tribution is labelled 214 Pb/ 214 Bi. Thorium, radium and 
lead isotopes can all be extracted with high efficien- 
cies from the underground water system by the HTiO 
ion-exchanger [10]. However, at the elution stage, very 
little thorium and lead are eluted by the 0.03M HN0 3 
solution. Thus, almost no 228 Th (which is above 224 Ra 
in the decay chain) can be found in the final counting 
samples, as illustrated in Fig. 4. 

During the time interval between the end of the elu- 
tion and the beginning of counting, which was nor- 
mally 25 hours for the hollow fibre procedure and 
5 hours for the pleated membrane procedure, the Pb 
and Bi isotopes will build up from the Ra isotopes. 
That is why contributions from the Pb and Bi isotopes 
are seen at the beginning of the counting period. (In 
the U chain their activities are relatively high in both 
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Fig. 4. The fitting of true ft — a coincidences binned per hour as 
the counting progressed in time for assay RDA5 for Th and U 
chain activities. The one sigma errors on the fitted parameters are 
statistical. 



background and assay spectra, because of 222 Rn con- 
tamination from exposure of the acidic eluate and liq- 
uid scintillator to the air.) Because of their short half- 
lives, Pb and Bi isotopes will decay away very fast and 
do not significantly affect the fitting accuracy for the 
Ra isotopes. The calculated activities (extracted plus 
background) of the relevant isotopes at the beginning 
of the counting are listed above Fig. 4. For the RDA5 
assay with hollow fibre filters, the background count- 
ing rates were 0.15 cph and 0.55 cph for the Th and 
U chain, respectively. 

Taking into account the decay of 224 Ra during the 
transport and processing time of the filters (for RDA5 
a factor of 0.78) and all of the relevant efficiencies 
(Table 1), the fit results for 224 Ra and 226 Ra can be 
converted to give the activity of each isotope on the 
filters at the end of extraction (given below in Table 3). 



4.2. Disequilibrium within the Thorium and 
Uranium chains 

Thorium tends to plate out on some surfaces but can 
form complexes in water, while radium can be leached 
out of surfaces in contact with the water. As a result 
it is possible that in an assay 228 Th (tjy 2 = l.9y) and 
its daughter 224 Ra (t^ 2 = 3.7d) are not in equilibrium 
by the time the water has reached the assay column. 
During the assay any thorium in the water will be 
extracted with ~95% extraction efficiency [10], and 
in the subsequent elution of the radium off the filters 
very little thorium is removed. So if the filters are left 
after being eluted after an assay, 224 Ra will grow in 
from the decay of any extracted 228 Th, and can then 
be eluted again and the activity measured. 

The result for the RDSA9 salt assay, and for an ear- 
lier assay (but with less precision) of pure heavy water 
in a storage tank, was that the Th and Ra activ- 
ities appeared to be in equilibrium, so that the 224 Ra 
would be supported during an assay. For the LWA5 
and LWA7 light water assays it was found that there 
was little (<15%) 228 Th activity in comparison with 
the 224 Ra activity extracted, suggesting that leaching 
of radium (or plating of thorium, or both) is more sig- 
nificant in the light water where the PMTs are im- 
mersed. For the results shown in Table 3, it is assumed 
that radium is supported in the heavy water assays but 
not in the light water assays. 

Allowing for any decay of 224 Ra during an assay 
gives the amount in the water provided that there is 
no significant plating of radium between the cavity 
and the assay column. Any plating is not expected 
to be significant with the water flow present in the 
polypropylene piping and has been measured to be 
less than 2%. Any contribution of 224 Ra from the pip- 
ing between the acrylic vessel and the assay filters 
has been estimated to be less than the equivalent of 
0.6xl0~ 15 g/g of 232 Th [12]. The final results are in- 
terpreted as the equivalent amounts of 232 Th and 238 U 
in heavy water or light water assuming that all the ra- 
dionuclides are in equilibrium. 

As the deuteron photodisintegration rate, which 
mimics the neutrino NC interaction, depends on the 
lead isotope activities (see Fig. 1), it is important to 
avoid any radon contamination both in the heavy wa- 
ter and in the light water. The N 2 cover gas together 



11 



Table 3 

The final results for some of the D 2 and H z O assays from October 1999 to November 2001; the one sigma errors given are statistical, 
the systematic errors are discussed in the text. 



Experiment 


Date 


Tonne 


T/m 2 


//min/m 2 


g Ti/m 2 /# 


Kr 15 g Th/g 


io-' 4 g U/g 


D 2 target 












<3.8 


<3.0 


RDA2 


13 Oct 99 


40 


8 


12 


0.5/H 


< 12 


0.28(1 ±0.76) 


RDA3 


03 Nov 99 


72 


15 


8 


0.5/H 




0.23(1 ±0.37) 


RDA4 


07 Dec 99 


209 


44 


16 


0.5/H 


5.0(1 ±0.14) 


0.03(1 ±0.42) 


RDA5 


22 Feb 00 


361 


75 


17 


0.5/H 


3.6(1 ±0.11) 


0.03(1 ±0.19) 


RDSA8 


30 May 01 


414 


100 


15 


2.2/P 


2.0(1 ±0.11) 


0.28(1 ±0.04) 


RDSA9 


24 Jul 01 


185 


44 


10 


2.5/P 


2.2(1 ±0.16) 


0.17(1 ±0.05) 


RDSA10 


25 Sep 01 


205 


49 


19 


2.5/P 


2.8(1 ±0.14) 


0.06(1 ±0.15) 


RDSA11 


23 Oct 01 


268 


64 


19 


2.5/P 


1.1(1 ±0.20) 


0.03(1 ±0.43) 


H 2 target 












<37 


<45 


LWA4 


16 May 01 


11.0 


21 


19 


2.5/P 


21(1 ±0.62) 


1.05(1 ±0.66) 


LWA5 


13 Jun 01 


18.6 


36 


19 


2.5/P 


79(1 ±0.09) 


0.48(1 ±0.25) 


LWA6 


26 Jun 01 


34.3 


66 


19 


2.5/P 


58(1 ±0.08) 


0.50(1 ±0.30) 


LWA7 


08 Sep 01 


27.5 


53 


19 


2.5/P 


85(1 ±0.06) 


0.31(1 ±0.17) 


LWA8 


06 Nov 01 


42.6 


82 


19 


2.5/P 


57(1 ±0.07) 


0.26(1 ±0.26) 



# H: Hollow fibre, P: Pleated membrane. 



with an adverse temperature gradient near the surface 
of both the D 2 and H 2 stops any significant radon 
contamination from mine air. The largest sources of 
radon in the D z O are from when the heavy water cir- 
culation system is turned on (bringing in radon from 
assay columns and from permeation through seals) 
and from when sources are deployed. 



5. Results and Conclusions 

The final results for some representative assays, af- 
ter subtraction of the measured background for each 
assay, together with the target limits for 232 Th and 238 U 
in the light and heavy water for SNO, are given in Ta- 
ble 3; also shown are the amount of water (tonne/m 2 ) 
and the flow (//min/m 2 ). It is assumed that there was 
no mixing of returned with sampled water in any of 
the assays. The activity is assumed to be uniformly 
distributed throughout the light and heavy water. For 
each of the assays with the pleated membrane filters, 
the measured backgrounds were similar and equivalent 
to a level of 8 x 1(T 16 g/g of 232 Th and of 5 x 1(T 16 g/g 
of 238 U for an assay of 200 tonnes of heavy water 

The one sigma errors given in Table 3 are statis- 
tical. The systematic error arises principally through 
the uncertainty of ~ 24% in the counting efficiency 



and the efficiencies of the chemical procedures (see 
Table 1), and to a lesser extent from the uncertain- 
ties in the volume of water sampled, in the amount 
of radium plating and in the background subtraction. 
The total systematic one sigma errors are estimated to 
be for the Th chain in D 2 +26—26% and in H 2 
+24-24%; for the U chain in D 2 +31-31% and 
in H 2 +30-30%. For the Th chain there is an addi- 
tional systematic error of 0.6x 10~ 15 g/g from the limit 
on any 224 Ra contribution from the piping. 

The light water assay LWA8 (see Table 2) shows 
that good extraction is achieved using the pleated 
membrane filters loaded with 2.5 g Ti/m 2 for ^80 
tonnes/m 2 and ~20 //min/m 2 . The assays LWA4-8 
suggest that the level of 232 Th is about 50 x 10~ 15 g 
per gram of light water. 

Assuming equilibrium in the decay chain and no 
plating of radium, it can be seen that the HTiO system 
has shown that the heavy water has met its limit of 
< 3.0 x 10~ 14 g 238 U per gram of heavy water and its 
limit of < 3.8 x 10~ 15 g 232 Th per gram of heavy wa- 
ter i.e. the Ra impurities in the heavy water are suffi- 
ciently low that the photodisintegration background to 
the predicted NC signal in the SNO experiment from 
the standard solar model flux of neutrinos is less than 
10%. However, while the light water has easily met 
its limit for the reduction of radioactivity backgrounds 
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of < 45 x lCT 14 g 238 U per gram of H z O it appears to 
be a little above its limit of < 37 x lCT 15 g 232 Th per 
gram of H 2 0. 

The total procedural efficiencies are 24±5% for the 
Th chain and 32±8% for the U chain. At present the 
total procedural backgrounds for both Th and U chains 
are low and at the level of 0.50 ±0.07 and 1.34 ±0.26 
counts per hour, respectively. The background errors 
include a contribution from the systematic error arising 
from the variation in the background rate from assay 
to assay. Defining the detection limit as 3 times the 
standard deviation of the background [19], the HTiO 
system currently has the sensitivity to measure down 
to - 3 x 10~ 16 g Th/g and to - 3 x 10~ 16 g U/g. 
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